Unitarity Bounds and the Cuspy Halo Problem 
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Conventional Cold Dark Matter cosmological models predict small scale structures, such as cuspy 
halos, which are in apparent conflict with observations. Several alternative scenarios based on 
modifying fundamental properties of the dark matter have been proposed. We show that general 
principles of quantum mechanics, in particular unitarity, imply interesting constraints on two pro- 
posals: coUisional dark matter proposed by Spergel & Steinhardt, and strongly annihilating dark 
matter proposed by Kaplinghat, Knox & Turner. Efficient scattering required in both implies m < 12 
GeV and m < 25 GeV respectively. The same arguments show that the strong annihilation in the 
second scenario implies the presence of significant elastic scattering, particularly for large enough 
masses. Recently, a variant of the coUisional scenario has been advocated to satisfy simultaneously 
constraints from dwarf galaxies to clusters, with a cross section that scales inversely with velocity. 
We show that this scenario likely involves super-elastic processes, and the associated kinetic en- 
ergy change must be taken into account when making predictions. Exceptions and implications for 
experimental searches are discussed. 
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I. INTRODUCTION 

There is a long history of efforts to constrain dark mat- 
ter properties from galactic structure (e.g. Recent 
numerical simulations [|| sharpen the predictions of Cold 
Dark Matter (CDM) structure formation models, and 
apparent discrepancies with the observed properties of 
structures from galactic to cluster scales are uncovered. 
The main one that has attracted a lot of attention is the 
cuspy halo problem , namely that CDM models predict 
halos that have a high density core or have an inner pro- 
file that is too steep compared to observations ( but 
see also Q). This has encouraged several proposals that 
dark matter might have properties different from those 
of conventional CDM (see g and summary therein). 

On the other hand, general principles of quantum me- 
chanics impose non-trivial constraints on some of these 
models. We focus here on the proposals of coUisional 
or strongly self-interacting dark matter (SIDM) by Jq] 
and of strongly annihilating dark matter (SADM) by [Q . 
Both require high level of interaction by particle physics 
standard: an elastic scattering cross-section of (Tel. ~ 
10"^'*(mx/ GeV) cm^ for the former and an annihilation 
cross-section of (Xann.Wrci ^ 10~^®(mx/ GeV)cm^ for the 
latter, where mx is the particle mass, and Vrei is the rel- 
ative velocity of approach. The proposed dark matter 
is therefore quite different from usual candidates such as 
the axion or neutralino. We show that the unitarity of 
the scattering matrix, together with a few reasonable as- 
sumptions, imposes interesting particle mass bounds as 
well as other physical constraints. This is done while 
making minimal assumptions about the nature of the 
interactions. Our results complement constraints from 
experiments or astrophysical considerations e.g. [B. 



Griest and Kamionkowski ||^ previously derived similar 
mass bounds related to the freeze-out density of thermal 
relics, assuming 2-body final states. In §||, we provide 
a general derivation for arbitrary final states using the 
classic optical theorem pO| ]. We summarize our findings 
in §111, discuss exceptions to our bounds, and other so- 



lutions to the cuspy halo problem. 



II. DERIVING THE UNITARITY BOUNDS 

Different versions of the unitarity bounds can be 
found in many textbooks, and can be most easily un- 
derstood using non-relativistic quantum mechanics (e.g. 
[prl ), which is probably adequate for our purpose. How- 
ever, the results and derivation given here might be of 
wider interest e.g. for estimating thermal relic density. 
Here we follow closely the field theory treatment of fl^ . 

The optical theorem is a powerful consequence of 
the unitarity of the scattering matrix S, i.e. S'^S ~ 1, 
which implies (1 - 5)t(l - 5) = (1 - S^) + (1 - S), or 



d7(/3|l-5|7)(7|l-^V> = 2Re(/3|l-5|a) (1) 



where a and /3 represent two specified states and 7 rep- 
resents a complete set of states with measure d'j. Using 
the definition of the scattering amplitude Af^a 

(/3|1 - S\a) = ^i{27t)^S\pp - p^)Ap^ (2) 

where pfs and pa are the total four-momenta, one obtains 



d7(2^)M^(pc. -p^)\A^c.V = 2ImA, 



(3) 
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if/3 = aineq. (||). We are interested in the case where a 
represents a 2-body state oi X + X ot X + X approach- 
ing each other. The final state 7, on the other hand, 
is completely general, and the integration over 7 covers 
the entire spectrum of possible final states. To be more 
precise, suppose \a) — |fci, sf ; A;2, Sj; n) where sf and 53 
represent the spin-states of the two incoming particles 
with spins si and S2 (in our particular case, si = S2) 
while fci and ^2 are their respective 4-momenta, and n 
labels the particle types (e.g. mass, etc.). Recalling that 
da/d'y oc |A^Q,p, eq. (^) gives in the center of mass frame 
(adopted hereafter i.e. ki -|- k2 = 0): 



^7^(« 



7) = 



Im Aaa 



2{Ei+E2)\ki 



(4) 



where the left hand side is exactly the total cross-section. 
This is the optical theorem. It states that the total cross- 
section for scattering from a two-body initial state to all 
possible final states equals the imaginary part of the two- 
body to two-body forward scattering amplitude. 

To use this theorem, we expand the scattering am- 
plitude in terms of partial waves i.e. states labeled as 
|ktot., £'tot., j, j^, ^, s, n) where ktot. is the total linear mo- 
mentum (= in the center of mass frame), -Etot. is the 
total energy , j is the total angular momentum, is its 
z-component, £ is the orbital momentum and s is the to- 
tal spin. Inserting appropriate complete sets of partial 
wave outer products into eq. (^, we obtain 

4i{27Tf[E[ + E!,/\k[\]i[Ei + E2/\ki\]^ (5) 

j.J^ i',s',£,s 

E {s's^'\sr4')s',..',{jfr'sni',A^i\'^'n 

where the crucial assumption is that S is rotationally 
invariant and so j and jz are conserved, in addition to 
energy conserving. The notation {t's'n'\\ — S\£sn)j_Etat 
emphasizes that S is diagonal in j, j'^ and i?tot but the 
dependence drops out because S commutes with JxiiJy 
The inner products (ss^|sf S2)si,s2 and {jj^\£^s^)i^s give 
the Clebsch- Gordon coefficients, and (ki|^£^) = ^^^-(ki) 
is the spherical harmonic function. We assume ki = z 
in which case y«.(ki) = (5f.,oV2^+ l/(47r). The index 
/3 denotes a 2-body final state \k'i, s'^; ^2,32; n'). 

Setting /3 = a, and averaging over the spin-states (i.e. 
(2si -I- l)-i(2s2 + l)"^Esj,.;) on both sides of eq. (|), 
the optical theorem, we obtain |l2| : 



2tt 



ftot. 



|ki|2(2si-f 1)(2S2 + 1) 

J 

E Re(^sn|l - S\lsn)j^E,^^ 



(6) 



This gives the total spin-averaged cross-section for scat- 
tering from X -I- X or X -I- X to all possible final states. 

For X + X annihilation, we exclude from the above the 
contribution due to elastic scattering (where type and 
mass of particles do not change i.e. X + X ^ X + X , 
implying |k'j^| = |ki|) To do so, we need the following 
expression for 2-body to 2-body scattering cross-section: 



da 



d(3 



4(Sl+S2)|kl| 



2n)^S\pf3-pa.)dP, 



(7) 



We average over initial spin states and integrate over out- 
going momenta, but focus on the elastic contribution (n' 
in 1/3) = \k[, s[; k'2, n') is set to n in jo;)) |12|: 



fel. 



|ki|2(2si + l)(2s2 + l) 

J 

E \{e's'n\l-S\£sn),^E.^, 



(8) 



The above is the total cross-section for elastic scattering 
(note: the same expression also describes X+X — » X+X 
elastic scattering) that has to be subtracted from txtot. to 
yield the total inelastic scattering cross-section, which is 
relevant for annihilation into all possible final states: 



CTincl. 



fc2(2si + l)(2s2 + l) 



(9) 
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1 - \{£sn\S\£s72)\'' - E \{e's'n\l~ S\£sn)\^' 



From eq. 



& (O), we can derive two bounds: 



^tot. <47r[|ki|2(2si + l)(2s2 + l)]-iEE2j + l (10) 

CTinel. <^[|kiP(2si + l)(2s2 + l)]-lEE2j + l (H) 

The first inequality uses |(^sn|S'|i?sn)p < 1, ob- 
tained from / d'y{£sn\S''\-f){-f\S\£sn) > \{£sn\S\£sn)\^ 
and S^S = 1. A similar bound can be derived for crd. 
as well, which coincides exactly with that for atot. ■ 

We pause to note that the above bounds assume only 
unitarity and the conservation of total energy and lin- 
ear and angular momentum. No assumptions are made 
about the nature of the particles, whether they are com- 
posite or point-like. Nor do we assume the number of 
particles in the final states. To obtain useful limits from 
the bounds, we take the low velocity limit. Assuming the 
scattering amplitude Apa is an analytic function of ki as 
ki — > (exceptions will be discussed in j |lll[ ) , and noting 
that k^{'k\££^) is a polynomial function of k, we expect 
the £ partial wave contribution to Af^a (eq. ^ to scale 
as |ki|^. This means in the low velocity limit, as is rele- 
vant for our purpose (typical velocity dispersion in halos 
range from 10 to 1000 km/s <C c), the ^ = or s-wave 
contribution dominates. Setting £ = in eq. (|lO|), (11): 
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0-tot. < l6Tr/{mxVrc\f 



CTinclWrol < 47r/(TO^Urel) (12) 



where = m\\v2 — vip/4 
second inequality agrees with 



™x^rci/4 is used. The 
|. Hence, 



(Jtot. < 1-76 X 10 



-17 2 

cm 



GeV 



mx 



lOkms 



-1 



O-inol.Wicl. < 1.5 X 10 



GeV 



mx 



lOkms-i 

Vrcl. 



(13) 
(14) 



Furthermore, if ainoi. is bounded from below, say a-mci. > 
Cann., One Can derive a lower bound on a^i. using eq. 
(§ & d), and setting £ = 0. Defining {X)j = [{2si + 
l)(2s2 + l)]"^ I]j,£,s(2j + 1)-'^, using S at the moment to 
denote {£sn\S\£sn), and noting that (1) j = 1 for £ = 0, it 
can be shown (7r/fc2)(l _ (|^|)2 ) > (7r/fc2) (l - (|5|2) j) > 

CTinei. > cTann., which impHcs (151) J < ^/T^^^a~J^^■ 
Also, del. > (Vfc?)(|l - S\^)j > (Vfc?)((l - \S\)')j > 
(7r/fc2)(l — (|5|),/)2. Combining, we have 



O-ol. > (7r/fcf)[l - Wl - fcfCTann./TI" 



(15) 



This tells us that the elastic scattering cross-section can- 
not be arbitrarily small given a non-vanishing inelastic 
cross-section, e.g. via annihilation. 

The above 3 bounds are the main results of this section. 
Two more results will be useful for our later discussions. 
For two-body to two-body processes, recall that the £, £' 
contribution to A^a scales as |ki|^|k'j^|^ . Using da/d^ = 
|A;3„|2(|ki|/|ki|)/[647r2(£;i -f £;2)2] (obtained from eq. @ 
by integrating over (3 except for solid angle fi), it can be 
seen that for clastic scattering, where |k'j^| = |ki|. 



da/dn const. [1 -I- 0(t;roi.)] 



(16) 



as |ki| 0. For inelastic scattering where the sys- 
tem gains kinetic energy by losing rest mass (e.g. de- 
excitation of a composite particle or annihilation), since 
jk'j^l approaches a non-zero value as |ki| — > 0, we have 



dcr/dfl — > ( const. /wi.oi.)[l + 0{vi-ci.)] 



(17) 



instead in the low velocity limit. The opposite case where 
the particle gains mass is discussed in |ll2| . 



III. DISCUSSION 

We can derive the following four constraints for 
strongly self-interacting dark matter (SIDM) |^ and 
strongly annihilating dark matter (SADM) ||^. 

1. The range fioi. '-^ 10^24 _ 10-23 cin2(mx/ GeV) is 
given by Q for SIDM to yield the desired halo properties. 
Using the lower CTol, and v^^i. 1000 km/s (appropriate 
for clusters), eq. (|f|) tells us mx < 12 GeV for SIDM. 

2. The annihilation cross-section from H, o'ann.'L'roi. ~ 
10~28 cm2(TOx/ GeV), together with eq. (p^) and v^d. ~ 



1000 km/s, gives us a bound of mx < 25 GeV for strongly 
annihilating dark matter. 

3. For SADM, efficient annihilation (a form of inelastic 
scattering) inevitably implies some elastic scattering as 
well. From eq. (p^), and using v^d ^ 1000 km/s as 
before, we have 



CTei. > 4 X lO-^^cm^lGeY/mxT 

[1 - Vl - 7 X 10-5(mx/GeV)3]2 



(18) 



Two simple limiting cases: when mx is close to the up- 
per bound of 25 GeV, a^. > 4 x 10^22 ^ni^; when mx 
is small, (Toi. > 5 x 10""^^ cm2(mj5s:/ GeV)**. Hence, elas- 
tic scattering is inevitable in this scenario, but can be 
reduced by having a sufficiently small mass. 

4. Recent simulations suggest that the simplest version 
of SIDM fails to match simultaneously the observed halo 
properties from dwarf galaxies to clusters ||l^,|| (see also 
|]l^ ), which have Uid. ranging over 3 orders of magnitude. 
It was suggested that an elastic scattering cross-section 
of (T oc 1/wroi. might solve the problem. But as shown 
in eq. (^), elastic scattering generally implies a — > con- 
stant in the small velocity limit. Hence, a oc l/wrd. likely 
requires inelastic processes. As eq. (17) shows, processes 
in which the net kinetic energy increases {\ki\ > |ki| 
in c.o.m. frame) can give such a velocity dependence. 
SADM provides an example. More generally, the net 
kinetic energy increase (super-elasticity) must be taken 
into account when considering the viability of a model 
with a cx l/t^rci. e.g. it may delay core collapse and make 
the core larger. Note, however, the general considerations 
in the last section does not forbid an elastic cross-section 
that increases as v,-ci. decreases e.g. the Oivid.) term in 
eq. ( |l6|) can have a negative coefficient. A 1/wrei. power- 
law may approximate such a cross-section, but likely only 
for a limited range of frci.. An example is the neutron- 
neutron scattering cross-section, which approaches a con- 
stant for |ki| < 10^2 GeV, and scales as l/wrci. only for 
10"^ < |ki| < 5 X 10-2 GeV 0. 

It is helpful to mention here possible exceptions to the 
above limits. Our bounds are obtained from eq. (|l3|) fc 
(|l^), which are the £ = (s-wave) versions of eq. ( [l0| ) 
& (|ll|). The argument for putting ^ = in the small 
velocity limit assumes the analyticity of at ki = 0. 
The latter breaks down if the interaction is long-ranged, 
e.g. Coulomb scattering. This is unlikely to be rele- 
vant, because there are strong constraints on dark matter 
with such long ranged interaction [Tsl . Our argument for 
the dominance of s-wave scattering can also be invalid if 
there is a resonance. However, given that the scattering 
cross section should vary smoothly over three orders of 
magnitude in velocities from dwarfs to clusters, a reso- 
nance seems unlikely. Finally, the most likely situation in 
which the bounds break down is if the particle has a large 
enough size, or the interaction has a large enough effec- 
tive range, R, such that |ki|i? > 1 (e.g. see [Esi). In such 
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cases, higher partial waves in addition to s-wave gener- 
ally contribute, and crtot. < 647ri?^ and our arguments 
turn into a limit on i? The condition |ki|i? > 1 gives 
the most stringent constraint on R for v-^d. = 10 km/s, as 
appropriate for dwarf galaxies: R > 10~^ cm(GeV/mjs: ). 
One can compare this with R for neutron-neutron scat- 
tering ~ 10"" cm 0. 

It is intriguing that halo structure might be telling us 
the elementary properties, in particular the mass, of dark 
matter. It is interesting that several proposals to ad- 
dress the cuspy halo problem, such as Warm Dark Matter 
p8[ and Fuzzy Dark Matter [|l^ make explicit assump- 
tions about the mass of the particles - mx ^ 1 keV and 
mx ~ 10"^^ eV respectively. For SIDM and SADM, as- 
trophysical considerations generally only put constraints 
on the cross-section per unit mass. We have shown here 
that unitarity arguments imply a rather modest mass for 
both scenarios as well. It is also worth pointing out that 
our arguments, with suitable modification to take into 
account bose enhancement and multiple incoming parti- 
cles, can be extended to cover dark matter in the form of 
a bose condensate, as has been proposed as yet another 
solution to the cuspy halo problem |2^]. They generally 
require small masses as well < 10 eV. 

A few issues are worth further investigation. Wan- 
delt et al. recently argued a version of SIDM, where 
the dark matter interacts strongly also with baryons, is 
experimentally viable, but requires mx > 10^ GeV, or 
mx < 0.5 GeV. Our bound here is inconsistent with 
the large mass region (but see exceptions above); exper- 
imental constraints on the low mass region will be very 
interesting (ctcI. < 10~^^ cm^). It would be useful to find 
a micro-physics realization of the coUisional scenario pl| ] 
or its variant where a scales appropriately with velocity 
to match observations. The impact of inelastic collisions 
on halo structures is worth exploring in more detail. It 
is also timely to reconsider possible astrophysical solu- 
tions to the cuspy halo problem, such as the use of mass 
loss mechanisms . We hope to examine some of these 
issues in the future. 
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